Melt viscosities of lattice polymers using a Kramers potential treatment. 
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Kramers relaxation times tk and relaxation times tr and tg for the end-to-end distances and 
for center of mass diffusion are calculated for dense systems of athermal lattice chains, tk is defined 
from the response of the radius of gyration to a Kramers potential which approximately describes 
the effect of a stationary shear flow. It is shown that within an intermediate range of chain lengths 
N the relaxation times tr and tk exhibit the same scaling with N, suggesting that iV-dependent 
melt- viscosities for non-entangled chains can be obtained from the Kramers equilibrium concept. 
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I. INTRODUCTION 

Extending the novel approach by KramersE! for poly- 
meric solutions under shear (based on replacing the sim- 
ple shear flow by its irrotational part)B~Q to melts al- 
lows us to compute a characteristic time tk related to 
the melt viscosity.Q In Ref. 5 the method was applied 
to alkane chains up to about 50 carbon atoms. In this 
note we compare this concept for melts with calculations 
of characteristic times for the end-to-end relaxation of 
chains (rotational time tr) and a time tq obtained from 
the center of mass diffusion. This is achieved by Monte 
Carlo simulations for dense systems of self-avoiding lat- 
tice chains on a simple-cubic lattice. Chains up to a 
length TV = 120 are considered, which near the upper 
limit already indicate the onset of entanglements. Our 
calculations show that there exists an intermediate range 
of chain lengths N where tk and tr agree in their depen- 
dence on N . This confirms the conclusion that important 
information on the melt viscosity as a function of N can 
be drawn from knowledge of fourth and sixth moments 
of the equilibrium distribution function for beads of the 
chain. In the next section we describe our lattice model, 
followed in section III by the description of our relaxation 
time results. 



II. LATTICE-MODEL AND KRAMERS 
RELAXATION TIME r K 

We model the self-avoiding polymer by athermal lat- 
tice chains where every point of a simple cubic lattice 
with unit spacing is allowed to be occupied by at most 
one chain monomer. In this way the hard-core repul- 
sion is taken into account. For the dynamics and for 
the equilibration of our system we use the generalized 
Verdier-algorithm,tj consisting of local moves of one or 
two monomers. This allows us to determine the center of 
mass coordinates, the center of mass diffusion constant 
and the end to end distance as a function of Monte-Carlo 
time steps, for a number of chain densities. The average 
chain densities are measured by the number of monomers 
per lattice point. We investigate monomer densities n 
ranging from a single chain up to n = 0.5 (half filled 
lattice). 



In Ref. 2 it was shown that a relaxation time tk can be 
obtained from the Kramers potential!!! for Rouse chains 
in a simple shear flow. This time is obtained from the 
equilibrium averaged fourth and sixth moments of chain 
monomer coordinates Xi and yi relative to the center of 
mass. If (Xq} is the mean squared x-component of the 
radius of gyration of the chain at equilibrium, thenQ 
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where Xi and yi are the coordinates of the chain 
monomers with respect to the center of mass of a chain. 
In the Rouse model one can calculate the shear viscosity 
•q from this relaxation time using 
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where p is the mass density, Na is Avogadro's number, M 
the molar mass of the molecule and r either the Rouse or 
the reduced Kramers time given by 
the reduced Kramers time tk, Eq. 
to give the melt (shear) viscosity of shorter alkanes up 
to about 50 CH2- units. In this note we do not obtain r\ 
directly but we calculate and compare tk , tr and tq as a 
function of chain length N at various monomer densities. 
While tk is calculated from (Q), we obtain tr by fitting 
an exponential to the decay of the end-to-end vector cor- 
relation function at large times. The time tq — (R G )/D 
is obtained from simulated center of mass diffusion coef- 
ficients D and the mean square radius of gyration (Rq)- 



III. SIMULATION RESULTS FOR tk, tr AND r G 

As a background to our further results we show in 
Fig. la plot of the diffusion constants D divided by the 
Rouse diffusion coefficient Do = (1 — n)/N for various 
densities. Here the factor 1 — n has been introduced to 
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account for the average effect of blocking of sites by other 
monomers. The drop in D/D a for sufficiently large A, 
seen in Fig. 1, can be thought to herald the onset of rep- 
tation. The horizontal line represents the expected A- 
dependences of the free chain while the line with slope 
— 1 reflects the expected behavior for reptating chains. 
The inset provides an approximate master curve for the 
diffusion constant versus A/A e where M e is an estimated 
length for the onset of entanglements B 

Fig. 2 presents double-logarithmic plots of the various 
relaxation times against A for three different densities. 
Naturally, for short chains one wouldn't expect any uni- 
versal properties so that differences in the behaviors of 
tk, Tr and tq are not surprising. In the dilute limit, 
displayed in Fig. 2a, both tk and tr grow with TV ap- 
proximately as [N — l) 2 - 26 when N > 20, ia-ceasonable 
agreement with earlier numerical studies. BEE!] For the 
purpose of illustrating this A-dependence common to tk 
and tr we have multiplied Tk by a constant factor so that 
both sets of data fall on a common curve. The associ- 
ated effective exponent in this range of A-values is in fair 
agreement with the theoretical expectation tr ~ 7V 2lAfl 
(see Ref. 11) with v w 0.59. For comparison, data for 
(R G )/N ~ N 2 "- 1 are also shown in Fig. 2a. On the 
other hand, with D ~ A -1 for unentangled chains, it is 
clear that tq will also be proportional to N 2v+1 . 

A different behavior is observed in Fig. 2b for a 
monomer density n = 0.3. Apart from a constant nu- 
merical factor, the relaxation times tk and tr fall to- 
gether in the range 10 < A < 10 2 . Deviations beyond 
A ~ 10 2 are expected to be caused by entanglement 
effects which then are no longer negligible and are not 
included in the Kramers method. Note that for large A 
sufficiently above the entanglement length one expects 
tr ~ A 3 4 and tq ~ A 3 . The effective dynamic expo- 
nent read from the slope of the straight line in Fig. 2b in 
principle may depend on n (see Ref. 10) but, similar to 
Fig. 2a, is still found in fair agreement with 2v + 1. On 
the other hand, screening effects now diminish the growth 
of (Rq) with A relative to the case of a free chain shown 
in Fig. 2a. The stronger increase of tq with A becomes 
clear from the behavior of D, see Fig. 1. 

These trends with n are continued in Fig. 2c where the 
range of agreement between tk and tr shrinks to about 
10 < A < 30. Fig. 2c also suggests that for A sub- 
stantially larger than A e tk, assumes again Rouse- like 
scaling behavior, i.e. tk ~ A 2 . 

In passing, we note that for the ratio of the square of 
the end-to-end distance R to the squared radius of gyra- 
tion Rq saturates with A and comes close to the value 6 
for our higher densities, as expected for Gaussian chains, 
whereas for lower densities this ratio is significantly larger 
than 6. This ratio versus A is plotted in Fig. 3. 



IV. CONCLUSION 

The Kramers potential approach for the viscosity re- 
laxation time of dense systems of lattice chains has been 
verified numerically for an intermediate range of chain 
lengths. For longer chains the breakdown of this concept 
could be identified as a result of the onset of entangle- 
ment effects. Generally, the Kramers approach offers the 
possibility to obtain information on a dynamic quantity 
from equilibrium data. Equilibrium lattice chain config- 
urations can be generated-,very efficiently via the con- 
figurational bias method. t2l From such data the static 
moments in Eq. 1 and hence information about melt vis- 
cosities in a certain range of densities and chain lengths 
can be obtained with minimal computational effort. 
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Fig. 1. Chain-length dependence of the diffusion coef- 
ficient for a free chain (□), n = 0.1 (*), n = 0.3 (x) and 
n=0.5 (+). By using Dq = k*/N with k* being chosen as 
a free parameter for each n and by rescaling the N axis 
via N — ► N/N e all curves (except the one for free chains) 
fall together on one curve, as shown in the inset. The val- 
ues chosen for N e are 23 for n — 0.5, 50 for n — 0.3 and 
150 for n = 0.1 (and k* = 0.50, 0.72, 1.28 respectively). 



Fig. 2. Comparison of the reduced Kramers relax- 
ation times tk/(N— l) 2 x with the relaxation time of 
the end-to-end vector r R /(N-l) 2 (□) and t g /(N-1) 2 
(*) for different number densities n, a) free chains, b) 
n = 0.3 and c) n = 0.5. Also included is the radius of gy- 
ration squared (+) divided by TV — 1. Typical error bars 
for larger iV— values are indicated in Fig. 2c, whereas the 
size of error bars in the intermediate range, where tr and 
tk coincide, is comparable to or smaller than the symbol 
size. 



Fig. 3. Comparison of R 2 /Rq vs. chain length N for 
different concentrations. Also included is the N — > oo 
limit for Gaussian chains. 
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